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Abstract 

We discuss SU(N) gluo-dynamics at finite temperature and on a spatial 
circle. We show that the effective action for the Polyakov Loop operator 
is a one dimensional gauged SU(N) principle chiral model with variables 
in the loop space and loop algebra of the gauge group. We find that the 
quantum states can be characterized by a discrete #-angle which appears 
with a particular 1-dimensional topological term in the effective action. We 
present an explicit computation of the partition function and obtain the 
spectrum of the model, together with its dependence on the discrete theta 
angle. We also present explicit formulae for 2-point correlators of Polyakov 
loop operators and an algorithm for computing all N-point correlators. 
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Yang-Mills theory in 1+1-dimensions is a prototypical example of a topo- 
logical field theory The classical theory is trivial and the quantum theory has 
no propagating degrees of freedom. It nevertheless has interesting features, 
particularly in the large N limit or when it is defined on a multiply connected 
spacetime [p], |] where its partition function depends on the area as well as 
topological invariants. On a sphere it exhibits a large N phase transition 
for some critical value of eA where e is the coupling and A is the area of 
the sphere ||. Wilson loops can be computed in both the strong and weak 
coupling phases and exhibit interesting behaviors [[§, g|. It has been shown 
that the strong coupling expansion can be rewritten as a lower-dimensional 
string theory ||. Furthermore, when defined on Riemann surfaces with 
non-zero genus it is known to have degrees of freedom related to the gauge 
group holonomy on the homology cycles of the surface 0], [j| . In particular, 
on cylindrical spacetime, it can be solved explicitly and has quantum me- 
chanical degrees of freedom corresponding to the eigenvalues of the Wilson 
loop operators which wind around the compact spatial direction 0-0 



In that case, topology also affects the gauge fixing, [pTT[| . Gauge invariant 
two dimensional models are also intimately connected with one dimensional 
integrable quantum systems [|1^, 

In this paper we shall examine Yang-Mills theory on a torus. It cor- 
responds to finite temperature 1+1- dimensional gluo-dynamics on a spatial 
circle. We show that, in any 1+1-dimensional gauge theory, correlators of 
Polyakov loop operators, which are the Wilson loop operators that wind 
around the periodic time, can be computed as correlators in a particular 1- 
dimensional gauged principle chiral model. In either pure Yang-Mills theory, 
or in Yang-Mills theory coupled to matter which is in the adjoint represen- 
tation of the gauge group, the effective field theory has a global symmetry 
which transforms the Polyakov loop operator by elements of the center of the 
gauge group. For the case of pure Yang-Mills theory, we find the effective 
field theory explicitly. We also find a solution of the effective theory and use 
it to give an explicit computation of the partition function and the two-point 
correlators of Polyakov loop operators as well as an algorithm for computing 
all correlators of Polyakov loop operators. 



A few years ago, Polyakov |L4| and Susskind [JT5J suggested that the 
question of confinement in quantum gluo-dynamics at finite temperature is 
intimately related to the realization of a global symmetry involving the center 
of the gauge group. This symmetry governs the expectation value of the 



Polyakov loop operator 

P(x) = trVexp U £ /T dTA (x,r)\ (1) 

which measures the gauge group holonomy in the periodic Matsubara time, 
<P{x 1 )...P{x m )P\y l )...P\y n ) >= 

_ J^e-f^^PW . ..P(x m )P\ yi ) . ..P\y n ) 



(The Euclidean path integral has periodic boundary conditions, A^l/T, x) = 
Afj,(0,x) and T is the temperature. For a discussion of the path integral 
formulation of finite-temperature gauge theory, see |16| .) The action and 
measure in the path integral (0) are symmetric under gauge transformations 

A M (r, x) -> g\r, x)A (t (r, x)g(r, x) + ig^r, a;)V^(r, x) (3) 

where g(r,x) is a gauge group element (here we consider SU(N)) which is 
periodic up to an element of the center of the gauge group, Z(N), 

g(l/T,x)=g(0,x)e 2mn / N (4) 

The path-ordered exponential transforms as 

tiVe l ^ /TdTAo ^ -> tig ] (l/T,x)Ve l ^ /TdTAo ^g(0,x) = ttPe* ^ drA ^ 

(5) 

The expectation value (Q) is interpreted as the free energy of the system in 
the presence of an array of external fundamental representation quark sources 
at positions x%, . . . , x m and anti-quarks at y\, . . . ,y n . The realization of the 
Z(N) symmetry is related to confinement - it is represented faithfully in a 
confined phase and it is spontaneously broken in a de-confined phase. In 3+1- 
dimensional chromodynamics a first order phase transition between confining 
and deconfining phases is expected at some critical temperature. The domain 
walls and Z(N) bubbles between different Zn vacua that should appear in the 



de-confined, broken Z{N) phase, have been extensively investigated [|17], [L£ 

On the other hand, 1+1-dimensional gluo-dynamics is confining and the 
global Z(N) symmetry is unbroken at all temperatures. There, the correla- 
tors of Polyakov loop operators are known to be related to the correlators in 



a generalized Sutherland model fT2|| . 



The Hamiltonian of 1+1 -dimensional gluo-dynamics is 



H= 6 4 [ L dxJ2E a (x) 2 (6) 
2 Jo V 

with nonvanishing commutator 

[A a (x),E b (y)]=i5 ab 5(x-y) (7) 

gauge constraint 

G a {x) = VE a {x) + f abc A\x)E c {x) ~ (8) 

and periodic boundary conditions, E a (L) = E a (0), A a (L) = A a (0). In 
Eqs.(||,^,|D a,b — 1 . . . N 2 — 1, the traceless hermitean generators T a are 
normalized so that trT a T b = \5 ab ) the structure constants are defined by 
[T a } T b ] = if abc T c and we define the N x N matrices E(x) = E a (x)T a } 
A{x) = A a (x)T a . 

Since A(x) and E(x) transform under the adjoint action of the gauge 
group, their periodic boundary conditions are preserved by time independent 
gauge transformations which are periodic up to an element of the center of 
the group, g n (L) = g n (0)e 2mn ^ N . G x = Jq dxG a (x)x a {x) commutes with the 
Hamiltonian and generates time independent, periodic gauge transforms 

e lG -E(x)e- iG - = gl(x)E(x)g (x) (9) 
e iG -A(x)e~ iG - = gl(x)A(x)g (x) + igl(x)Vg (x) (10) 

where go(x) = 1 — ix{x) + • • • and where go(0) = go(L). The constraint 
d§D indicates that the physical states are invariant under all strictly peri- 
odic gauge transforms. The coset group of all gauge transformations modulo 
periodic ones is isomorphic to the center of the gauge group, Z^. The phys- 
ical states must transform under an irreducible unitary representation of the 
coset. Representations of are one dimensional and are characterized by 
a discrete angle.f] 

2 The number of connected components of the group of periodic gauge transformations, 
which is the loop group based on SU(N)/Zn, is 

IT (loops on SU(N)/Z N ) = U^SV '(N) / 'Z N ) = T1 (Z N ) = Z N (11) 

The group Zn has N 1-dimensional irreducible representations, labelled by an integer 
m = 0, 1, . . . , N — 1 and where the elements are U (n; 9 m ) = e lSmn and 6 m = 2irm/N . 



This theta-angle will also characterize the quantization of Yang-Mills the- 
ory with coupling to any adjoint representation matter fields. Some related 
issues have been pointed out in jni. Also, for adjoint QCD on a hyper- 



torus in d spatial dimensions, there are d angles related to periodicity of 
gauge transformations as well as perhaps others related to higher dimen- 
sional winding numbers. In either case, if one introduces matter fields which 
are in the fundamental representation, their boundary conditions reduce the 
gauge symmetry to gauge transformations with a fixed periodicity, i.e. one 
of the degenerate states, and there is no theta-angle. 

In the Schrodinger picture, where states are wavefunctionals ^[A], and 
the electric field is realized as E a (x) = | SA t^ , the physical state condition 
is 

frmZ + ™ x) mz) »*-M = ° (12) 

The wavefunctionals of physical states therefore transform as 

% hys [giAg n + ig]yg n] 9] = e in0 *[A; 9} (13) 

The thermodynamic partition function is obtained by taking a trace of the 
Botlzmann factor, e~ H ^ T , over the physical states. This is most conveniently 
expressed as a trace using eigenstates of A(x) and a projection onto gauge 
invariant theta-states, 

Z = VOL~G ^ Ji d 9n]e~ idn jj[dA\x) <A\e- H ' T W n Ag n + ig]yg n > 

(14) 

The integration over gauge transforms projects onto physical states with 
angle 9. This integral can be expressed in terms of the standard Euclidean 
path integral (see |17| ) which is used to take expectation values of observables 
as in (0) by first considering the integral with twisted boundary conditions, 

z = votg £ / [^]e- ien n^ a (^^)^ /ol/Tdrd ^ trA2(x ' T) (is) 

A{l/T, x) = gl(x)A(0, x)g n {x) + ig\{x)Vg n {x) (16) 

The boundary conditions can be made periodic by performing a non-periodic 
gauge transformation. 

A(t, x) -> gl(r, x)A(t, x)g n (r, x) + igl(r, x)Vg n (T, x) (17) 



where g n (l/T,x) = g n (0,x)g n (x). The result is the standard form of the 
path integral in (H). The loop operator is identified as 



and the trace of the Polyakov loop operator is equal to the trace of the gauge 
group valued variable in (JTJ). 



P{x) = trg n (x) (19) 



The integrand in (0), 



e - id n < A \ e -H/T WnAgn + iglygn >s e -s es[ A, 9n] ±_ tl ( gn{ L)gi(O)Y 0N/27T 

(20) 

is a functional of A and g n which has the symmetry properties 

S e g[A, g n ] = S e $[v}Au + iu^Vu, vT l g n u] (21) 

S cS [A, g n z] = S eS [A, g n ] (22) 

The variables A and g n take values in the loop algebra and loop group of the 
gauge group, respectively. This is a 1-dimensional gauged principle chiral 
model. The theta-term in (|20|) is a topological term in the sense that it is 
insensitive to local variations of g n (x) and depends only on its holonomy on 
the circle. (p2|) implies that the effective action is invariant under the Zn 
symmetry g n — > g n z. The expectation values of Polyakov loop operators are 
given by the correlators of traces of the group elements 



<P(x 1 )...P t ( 2/fc ) >= 

-8N/2TT ( 23 ) 



ZnfdA(x)[dg n {x)]e- s ^ A >^Ui (g n (L)gl(oj) ^ txg n { Xl ) . . .txgl{y k ) 



EnIdA(x)[dg n (x)]e- s ^ A ^^tr (g n (L)gt(0) 



The equation (^) defines an d-dimensional effective action for the Polyakov 
loop operator in any d+ 1-dimensional gauge theory coupled to matter fields 
where there is a trace over the matter degrees of freedom on the right hand 
side of the formula. In all cases, the resulting effective action has the gauge 
invariance property in ([H]). When the matter is in the adjoint representa- 
tion, or any other representation which is invariant under transformations 



in the center of the gauge group, the effective action also has the global 
symmetry in d22|) . In the present paper, we shall concentrate on finding the 
effective action for pure 1+1-dimensional Yang-Mills theory. 
By evaluating the matrix element in fllj 



< A\ exp I— £ J Q dxjj^ j \glAg n + ig\ Vg n > (24) 



( T r L 

exp ( -— J dx tr(A - glAg n - ig ] n Vg r , 



(25) 



(when defined using zeta function regularization, the constant factor which 
should appear on the right-hand-side is one) we present the partition function 

as 

Z[9,T,L] = — / Y[dA a (x)[dg n (x)} 

exp (-^ £dx tr\Dg n \ 2 j tr (g n (L)gl (0))"^ (26) 

where Dg n = Vg n + [A, g n \. 

The integral in ( p6|) can be done by using gauge symmetry to diagonalize 
the unitary matrix g n = Ug^U' 1 where g^ix) = diagona^e 4 ^™^, e l ^ x \ . . . , e l< ^ n ( x >) 
The variables are angles G [0, 27r]. The measure in the integral has the 
form 

[d 9n (x)]= n n€wnj^(cw-^(^))cc)[^)] ,w) 

xG[0,L] a a <l3 Z a 

and the action is J / L feu W°V< + E a ,/3 | A^| 2 |e iC - | 2 ) . The in- 
tegral over A a p, where a ^ (3, cancels the Vandermonde determinant in the 
integration measure. The integration over diagonal components of A yields 
an infinite factor which compensates the infinite normalization of the plane- 
wave state which was used to take the trace in (14). The partition function 
is an integral over the eigenvalue variables, 



N-1 , N-1 



z = Ticvr E e ~ t9n I II II <^0*0 ex P " T / e2 / dx E VC VC + e wj 



N-1 / N-1 



(28) 



where the integration variables have the boundary condition <f)%(L) = 0°(O) + 
On + 2iTp a , 9 = 2-kitl/N and p G ZZ and the integration measure integrates 
4>(x), at each point x, over the range [0, 2ir) and sums over the integers p a . 
The integral is invariant under the field translation symmetry 



(29) 



where c a is a constant. Note that this is not a symmetry of the action or the 
integration measure in the path integral fl26| ) separately but only appears 
after the integration over A. This symmetry can be used to extend the 
limits on the integration over <p(x) at each x to infinity. One can untwist 
the boundary conditions and take into account the periodicity of the 
variables, by changing the functional integration variables to <f) a (x), by means 
of 

2n(n + Np )x 



(*)= v NL ' +r(x) , 

is real and periodic, (t> a {L) 



(30) 



with p a £ Z. In (|30|) (f) a (x) is real and periodic, (j) a (L) = (f) a (0). The parti- 
tion function can be written as the product of an analytic part Z a containing 
the integral over the periodic fields 4>(x) and a topological part Z t containing 



the summation over integers p a , Z = Z a Z t . Z a 
of a zeta function regularization. Expanding C 



can be computed by means 
(x) in modes 



r(x) 



V J-i k=-oo 



i2irkx/L 



(31) 



the functional measure is defined as Y\ x d(j){x) = T[ k dak- The integration 
over the zero mode produces an infinite, irrelevant, temperature independent 
factor proportional to the volume of the gauge group. The functional integral 
in Z a is proportional to the determinant of the laplacian, 



Za = JJ ^det ^W) 2 = det m e^-iX'W 



where Q is the (N — 1) x (AT — 1) matrix 



2T 



n 



4vrT 



( 2 


1 . 


• M 


1 


2 . 


. i 


Vi 


1 . 


. 2) 



(32) 



(33) 



and we have used zeta-function regularization with 

1 



OO 1 

k=X 



(34) 



and C(0) = -1/2 and C'(0) = -(log2?r)/2. 

Z t is given by a finite sum over n and an infinite sum over all the p a . 
Rewriting 



cxp 



i2nmn\ 



N 



1 



exp 



N-l 



i2nm ^ (n + Np a )/N 



a=l 



(35) 



the finite sum over n and the N — 1 infinite sums over the p a become just 
N — 1 infinite sums over the integers n a = n + iVp a . The ra a in fact span Zj 
when < n < N — 1 and p e Zj. Z t then reads 



Z t — ^ exp yZirin 1 ■ z + mftfQnj = 6 (z, Q) 



(36) 



Here we introduced the ^-function of several variables according to the con- 
ventions of Ref. p0| , ft, z are thought of as iV — 1-dimensional column vectors, 
n* = (ni, n 2 , • ■ • , njv-i), = (m/N, m/N, ■ ■ ■ , m/N). 

In order to obtain the energy spectrum of the theory we can now apply 
the generalized Poisson resummation formula 



exp y—ix gijrri m? — 2irim' l a, i 

rh=— oo 

where g JJ is the inverse of gij 



1 



det gij n 



exp -7rg lj (rii- ai)(nj - dj) 

=— oo 

(37) 



e 2 LN 
4ttT 



/iV-1 -1 
-1 JV-1 



V 



- 1 \ 
-l 

iV- 1/ 



(38) 



and det g^ = (4irT/e 2 lN 2 y N ^N. The Jdetgij cancels the temperature 
dependence of Z a so that the partition function reads 

v 21 



+oo 
n=— oo 



4T 



AT-l 



£ (iVn a + 



m — 



a=l 



iV 



■ (39) 



From ( |5SD one can deduce the energy spectrum of the theory for each m- 
sector. The result is energy levels which are characterized by iV — 1 integers 
n and 



E(n) 



4 



JV-l 



E ( Nn * + 



m) 



a=l 



1 

N 



rN-i 



E (Nn a + 



m) 



(40) 



which agrees with the one obtained, in a different context, in Refs. |7|, ||, |9| 
with the small difference that one finds that each energy eigenvalue belongs 
to a specific topological sector, labelled by m. For example for SU(2) in 
Refs. 0, H H was found £7 n = e 2 Ln 2 /S and here we see that the eigenvalues 
corresponding to even integers belong to the m = sector (E n = e 2 L(2n) 2 /8) 
those corresponding to odd integers (E n = e 2 L(2n + l) 2 /8) belong to the 
m = 1 sector. 

The calculation of correlators of Polyakov loop operators straightforward. 
The path integral is easily performed and the sums on n and the p a can be 
treated as before. 

Correlators are nonvanishing only if they are invariant under the trans- 
formations by the center of the group, Z(N), in d22|). This occurs when the 
number of g's in a correlator is equal to the number of g^s, modulo N. Using 
Eqs.(p7[), the 2-point correlator Pffl(x,y) =< tr{g(x)}tr{g^(y)} > m /N 2 can 
be presented as 



1 



ZN 2 



+oo 

E 



exp 



N-l 



+ E ex p 

where 
For SU(2) 



-■ng' 



m 

rij H 

1 N 



-7ig lJ 

x-y 

NL 



m 



TIa H 

3 N 



x-y 
NL 

m x — y 



m 

U3 + N 



NL 5ja 



G(x, y) = exp 



Le 2 (N - 1) (\x- y\ (x - yf 



ATN 



L 2 



x-y 
NL 

G(x,y) ,(41) 
(42) 



2Z 



E ex p 



e 2 L 



8T 



(2n + m) 2 - (An + 2m) (— 



\x - y\ 
L 
(43) 

Any other correlator of Polyakov loop operators can be analogously obtained. 
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